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Abstract 

In this paper we consider Wannier functions of quasi-periodic g-gap 
{g ^ 1) potentials and investigate their main properties. In particular, 
we discuss the problem of averaging underlying the definition of Wan- 
nier functions for both periodic and quasi-periodic potentials and express 
Bloch functions and quasi-momenta in terms of hyperelliptic a functions. 
Using this approach we derive a power series expansion of the Wannier 
function for quasi-periodic potentials valid at |a::| ~ and an asymptotic 
expansion valid at large distance. These functions are important for a 
number of applied problems. 

1 Introduction 

The Schrodinger operator with periodic potentials is characterized by two dis- 
tinguished complete sets of orthogonal functions: the Floquet-Bloch functions, 
introduced by Floquet (1882) and F. Bloch (1928), and Wannier functions, in- 
troduced by G. Wannier (1937). The Schrodinger operator with algebraically 
integrable potentials allows one to construct Bloch ( Floquet-Bloch in one di- 
mension) and Wannier functions and describe their properties with unprece- 
dented completeness. While Bloch functions were studied for a long time and 
there are many reviews and monographs on the topic, the study of the Wannier 
functions corresponding to finite-gap potentials has begun only recently |1(J| . 

The aim of the paper is to further expand these studies by presenting a 
definition and a complete description of the properties of Wannier functions for 
g-gap potentials with g > 1. To this regard, we remark that the usual Wan- 
nier functions considered in Condensed Matter Physics are defined for periodic 
potentials. By considering Wannier functions for (?-gap potentials, which are 
in general quasi-periodic functions, we are generalizing the usual Wannier func- 
tions to quasi crystals i.e. to a quasi periodic lattices, this being a first main 

'Institute of Magnetism, Vernadski str. 36, Kiev-142, Ukraine. 

^Department of Mathematics and Statistics, Concordia University, 7141 Sherbrooke West, 
Montreal H4B 1R6, Quebec, Canada. 

t Dipartimento di Fisica c Istituto Nazionale di Fisica della Materia (INFM), via S.AUende, 
184081 Baronissi (SA), Italy. 



1 



result of the paper. This circumstance leads to another interesting problem. 
As well known, in order to define Wannier functions one must use an averag- 
ing procedure which differs essentially for the periodic and quasi-periodic cases. 
Although the averaging problem has been considered for some algebraically in- 
tegrable system (see e.g. [21], ES], [33): to the best of our knowledge, it has 
not been discussed with sufficient completeness and not linked to Wannier func- 
tions, this being indeed a second main result of the paper. The study of the 
averaging problem in connection with integrable systems, however, goes beyond 
the frames of this paper. Here we shall only remind that averaging procedures 
also arise in Krylov-Bogoliubov-Mitropolsky perturbation theory of integrable 
systems, in the study of Seiberg-Witten theory, in multi-matrix models etc. We 
hope to expand this direction of our study in the future. 

The theory of hypcrcUiptic curves of different genera leads to a rich struc- 
ture of objects which inherits the rich structure of the moduli spaces of the 
corresponding hyperelliptic curves. For example, in the case of elliptic curves 
g = 1, periodic and ergodic cases coincide because there exists only one fre- 
quency. In the case of the two-gap potential the spectral variety is given by the 
ultra-elliptic g = 2 curve and periodic and ergodic cases are rather different. 
Moreover, periodic case admits further specialization to the elliptic periodic- 
ity. In the cases <? > 3 a new phenomenon appears: the curve admits so-called 
singular half-periods which leads to certain complications of the theory. 

Technically speaking, our development is based on the realization of hyper- 
elliptic functions in terms of multi-variable tr-function. This realization repre- 
sents a natural generalization of the Weierstrass theory of elliptic functions to 
hyperelliptic functions of higher genera. Higher genera theory was developed by 
K. Weierstrass and F.Klein and its exposition is fixed in classical monograph of 
H.Baker [2]; recent results in the area are given in |12II13[[T^I15| whilst various 
applications see e.g. in |TH ] . pS ] . P ] .|51.[H|. 

The restricted length of the paper does not permit us to give proofs of the 
theorems for which we plan a separate publication. 

The paper is organized as follows. In the Section 2 we describe quasi-periodic 
potentials from the pure spectral point of view. We show that the basic objects 
like finite-gap potential, Weyl and Bloch functions, which are defined initially 
on the complex plane can be reasonably lift to hyperelliptic Riemann surfaces 
whose genus coincides with the number of gaps in the spectrum. In the Section 
3 we collect necessary results from the theory of hyperelliptic cr-function such as 
hyperelliptic p and ^-functions and differential relations between then, as well 
as, addition theorems. We also describe certain sub-varieties of the 0-divisor. 
In Section 4 we show how to evaluate Bloch function and quasi-momentum in 
terms of multidimensional a functions. In section 5 we discuss periodic and 
quasi-periodic potentials and show how to compute averages of squared Bloch 
functions in both cases. In Section 6 we discuss Wannier functions for periodic 
and quasi-periodic potentials and derive a power series expansion of the Wan- 
nier function for quasi-periodic potentials valid at |a;| ~ and an asymptotic 
expansion valid at large distance. Finally in Section 7 the main conclusion and 
perspective of future work are presented. 
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2 The Schrodinger operator with quasi-periodic 
and finite-gap potential 

2.1 Quasi-periodic potential 

The function u{x) from a Banach space Cfc(]R) of continuous bounded functions 
is called almost periodic if the set {Txu{-), x G M}, where Txu{-) — u{- + x), is 
relatively compact in Cf,(M). A closure F of this set is a compact in metrizable 
Abelian group. On the set F there exists a normalized Haar measure /i which 
turns out to be T^-invariant and ergodic. Thus, each almost-periodic function 
generates a probability space (F, ^,Tj.). The operation of averaging on this 
space is given by 

ifiu)) = \im- f f{T,u)dx = f f{u)ti{du). (2.1) 



By means of the differential expression 

L{u) = -dxx + u{x) 

we define for w g F a Schrodinger operator in £<^(R), which is essentially self- 
adjoint. Let A € C and c(x, A), s(a;, A) denote solutions of the equation Llp = Xip 
with the initial data c(0, A) = 1, c'(0,A) = 0, s(0, A) = 0, s'(0,A) = 1. The 
functions c{x,X), s{x,X) are integral functions of order 1/2 with respect to A 
and continuous with respect to all variables (cc, A). The limits 

w±{X) ^ w±{X,u{-)) = T lim 4^ (2.2) 

x^oo S(X, X) 

exist and are called the Weyl functions. Besides these, we shall also employ 
the Weyl functions ■w±{x,X) = w±{X,Txu{-)). It is well known that all Weyl 
functions are holomorphic in C+ = {A S C | ImA > 0}, map C+ — *■ C+, and if 
they have zeros or poles on the real axis, these could be simple only |16| . By 
means of the Weyl functions wc can define the functions 

^p{x, X) = c{x, X) ± w±{x, X)s{x, X) = cxp ^± J w±{y,X)dy^, (2-3) 

which, for each A G C+, belong to iL^(M±), where IR.+ ~ [0, oo), and M_ = 
(— cxD, 0), respectively. By definition, the functions ^'±(2;, A) satisfy the equation 

-a,, V± (x, A) + (uix) - A)^± {x,X)= 0. (2.4) 

By substituting (|2.3f) into this equation, we have the following equation for the 
Weyl functions: 

±dxw±{x, X) + wl{X,x) + X - u{x) = 0. (2.5) 



3 



Following I2H1) we define the Floquet function 

f{x,X) = ^{w+ix,X)+w-{x,X)). (2.6) 
Let A = ^ + i?7. The Floquet function has a finite limit almost everywhere at 

^ e K, 77 i 

/(e + iO,x) =-/(x,e)+i7rn(x,0, (2.7) 

where is the Lyapunov exponent, and n(^) is the number of states. 

The number of states, n(a;,^), determines the spectrum S(?x) of the operator 
L(u): for a.e. u G F the spectrum is a set of growth points of the number of 
states PHI) 

Ti{u) = supp(c?n). 

The Lyapunov exponent ^(a;,^) determines an absolutely continuous spectrum 
i.e. a closed set of spectral points with zero Lyapunov exponent, 

^aAu) = {eeK|/(x,C) = 0}, 

where the closure is with respect to Lebesgue measure ESj • We now proceed 
to discuss a special subset of quasi-periodic potentials, the so called finite-gap 
potentials (more details see in 02] , > ESI ) ■ 

2.2 The finite-gap potential 

Theory of finite-gap potentials of the Schrodinger operator has long history 
which goes up to Hermite, Halphene, Darboux. Recent development was stim- 
ulated by the soliton theory and achieved in the works made at the middle of 
seventieth by Novikov, Dubrovin and Novikov, Its and Matveev and Krichever 
(see e.g. 021 and references therein). 

Definition 2.1 The almost-periodic function u{x) is called a finite-gap poten- 
tial if the spectrum of the Schrodinger operator L{u) ~ —dxx + u{x) is the union 
of a finite set of segments of Lebesgue (double absolutely continuous) spectrum 
£2] U [£^3, E4] U • • • U [-B231 , 00], where boundaries of the bands are supposed 
to be real and ordered as: Ei < E2 . ■ ■ , < £'2g+i < +00. 

Starting from this definition, we can derive explicit expressions and basic prop- 
erties of finite-gap potentials, Weyl and Bloch functions. We shall show that 
the adequate language to describe these object is the theory of hyperelliptic 
Riemann surfaces. To this end we fix the Riemann surface X of genus g of the 
algebraic curve given by the equation 

23-1-1 

fj.^ = R2g+i{X), i?23+i =4 [| (A-Sfe). (2.8) 

fe=i 

We denote the running coordinate of the curve X a.s P = (A, /i) and coordinate 
of branch points as {Ek, 0), k — 1, . . . , 2g + 1. 



4 



Theorem 2.1 The Weyl function of the finite- gap potential u{x) is defined on 
the Riemann surface X of the hyperelliptic curve 

= i?2s+i(A), i?23+i(A) =4 [] (A-S,) (2.9) 

i=l 

hy the formula 
where 

a 

5(a;,A) = n(A-Afc(a;)) (2.11) 
fe=i 

and points (Afe(x), ,Ufe(x)), k = 1, . . . ,g, are distinct points of the curve X eval- 
uated accordingly with changing of x. 

For our further discussion we need the expression of the number of states, n(^), 
in terms of the function Wj^{x, ^), 

ntt) = I = I (to »+(x,«)) = lf (i;;:^) <«. 

For a finite-gap potential u{x) we have 
Therefore in this case the number of states is 



The number of states and the wave number for all finite-gap potentials are 
connected as follows 

n(o = -m). 

A direct consequence of the above discussion is the following theorem. 

Theorem 2.2 The quasi-momentum of eigenfunction can be defined by the ex- 
pression 

*(p) = «.,p)) = (^)^/M,., p,.) 

where P = (A, /x) and Pq = (Ao,/io) o-re points of the hyperelliptic curve fj? = 
-R2a+i(A). 
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The last integral in Eq. 12.12|l is the Schwarz-Cristoffel integral which maps 
the l.h.s. half-plane of complex A-plane, to a rectangle in the complex /c-plane 
123 EHl- We shall further interpret this integral as a second kind Abelian integral 
on X over meromorphic differential, 

dOC = ^^^^l^^dA, (2.13) 
M 

which we shall call differential of quasi-momentum. 

Theorem 2.3 The eigenfunction of the g—gap potential, which is normalized 
by the condition 

(|V(x,P)p) = l, (2.14) 

is of the form 



{S{x,P)) 2 7 S{y,P) 

We see that in order to define the quasi-momentum and normalize the eigenfunc- 
tion we must fulfill the averaging. As a consequence we should say some words 
about the averaging for quasi- and almost-periodic functions since the potentials 
with finite number of gaps are quasi-periodic functions and the potentials with 
infinite number of gaps are almost-periodic functions. 

2.3 Averaging for the finite-gap potentials: periodic and 
ergodic cases 

Any almost-periodic function has a mean value 

1 



(/> - lim - / fix)dx. 
L^oo L Jq 

This allows for any almost-periodic function to build a Fourier series 
f{x) ~ ^ A„ exp(iA„x), An = (/(x) exp(iA„a;)). 

n 

The numbers A„ are designated as the Fourier frequencies and the numbers An 
as the Fourier coefficients of the function f{x). 

We remind that a countable set of real numbers Ai, . . . , A„, . . . has a ratio- 
nal basis ai,a2, ■ ■ ■ if the numbers ai, a2, . . . are linearly independent and any 
number A„ can be presented as their finite linear combination with rational 
coefficients, i.e. 

s„ 

A„ = ^r^a,, r^eQ. 

k=l 
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Wc say that the basis is finite if it is finite set, we say that the basis is integer 
if all numbers are integer numbers. If the Fourier frequencies of almost- 
periodic function have a finite and integer basis we designate the appropriate 
almost-periodic function as the quasi-periodic one. A quasi-periodic function 
with one period is pure periodic one. 

Function F{xi, x^, ■ ■ ■) of finite or countable set of variables, each of which 
admits all real values, is called limiting periodic if it is a uniform limit of periodic 
ones, i.e. if for any positive real number e we can point out such an integer 
positive number n{e) and such a periodic function Fs{xi,X2, ■ ■ ■ ,a;„) that 

sup \F{xi,X2, ■■■) - Fs{xi,X2, ■ . ■,Xn)\ < £■ 

— oo<Xfc,fc=l,2,...<-|-oo 

It appears that for any almost-periodic function /(x) there exists such a limiting 
periodic function F{x\,X2, ■ ■ ■) of finite or countable set of variables such that 

f{x) = F{x,x,...). 

In other words any almost-periodic function is a diagonal restriction of some 
limiting periodic function. The properties of the limiting periodic function 
F(a;i,a;2, . . .) depend essentially on the basis of the Fourier frequencies of the 
function /(.x). If the basis ai, a2. ... of the almost-periodic function /(.x) is in- 
teger then the limiting periodic function F{xi,X2, ■ ■ ■) is periodic with periods 

27r/ai, 27r/a2, If the basis a\,a2, ... of the almost-periodic function f{x) 

is finite then the limiting periodic fimction F{xi, X2, ■ ■ ■) depends on finite set 
of variables. If the basis ai, 0:2, • • • of the almost-periodic function f(x) is fi- 
nite and integer then the limiting periodic function F{x\,X2, ■ ■ ■) is a periodic 
function on finite set of variables. Obviously, the mean value of the function 
f{x) depends essentially on whether this function is periodic, quasi-periodic or 
almost-periodic. If the function f{x) is periodic then the mean value of the func- 
tion is just the average over the period. If the function f{x) is quasi-periodic 
or almost-periodic then the mean value of the function is the average over the 
phase space, i.e. ergodicity takes place. This is illustrated with the following 
classical statement. 

Theorem 2.4 (Kronecker— Weyl) Let Afe, k = 1, . . . ,n be real linearly inde- 
pendent numbers, 9k, k = 1, . . . ,n be arbitrary real numbers, dk,k = 1, . . . ,n be 
arbitrary positive numbers. Let x{xi,X2, ■ ■ ■ ,Xn) be a characteristic function of 
parallelepiped in K" defined by inequalities 

0k-h<Xk<dk + 6k, k = l,...,n. 

Then, uniformly in L we have 

lim / xi^ix-0i,...,XnX-9n)dx = -n:~'^di...5n, 

where the function x(xi, X2, . . . , a;„) is the periodic continuation to the whole 
M" in all variables Xk,k = 1, . . . ,n, with periods 27r. 
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In the problem under consideration a finite-gap potential appears as a re- 
striction to a straight line of some meromorphic function defined on a Jacobian. 
Of course on this line the potential must be real and non-singular. Depending 
on a direction of the line with respect to periods of the Jacobian this finite- 
gap potential can be periodic or quasi-periodic and this circumstance influ- 
ences the averaging essentially. The averaging of quasi-periodic functions which 
are restrictions of the 6'— functions was previously described in the literature 
(see e.g. 01]), and the averaging of solutions of integrable evolution equations 
(see e.g.|2I],|2ni)- We shall use these results when discussing the averaging for 
the finite-gap potentials. We first shall formulate our problem in terms of the 
tr— functions. 

We remark that it is convenient to use the language of cr— function in our 
foregoing discussions since all expressions in cr— functions are form-invariant with 
respect to the number of gaps of potential (or the genus of curve) and this fact 
is very useful for exposition. Keeping this in mind, we give in the next section 
a brief account of necessary facts from the theory of hyperelliptic ct— functions. 



3 Hyperelliptic cr— functions 

The theory of hyperelliptic cr— functions represents a many-dimensional gener- 
alization of the theory of elliptic functions and provide suitable language of our 
development. 

3.1 Curve and differentials 

Let X be the hyperelliptic curve given by the equation 

2g 2g+l 

^2 ^ 4^2g+l ^ ^ = 4 [J (A - i?fe) - R{\), (3.1) 

1=0 k=l 

realized as a two sheeted covering over the Riemann sphere branched in the real 
points (£^fc, 0), /c e S = {1, . . . , 2^ -I- 1}, with Ej ^ Ek for j ^ k, and at infinity, 
E'2g+2 = oo. The order of branch points is according to Ei < E2 < ■ ■ ■ < E2g+i 
(see Fig.lH). 

A set of 2g canonic holomorphic differentials dh = (d/ii, . . . , d/ig)^ and 
meromorphic differentials dr — (dri, . . . , drg)^ is defined by the following ex- 
pressions, 

dh, = dA, j = l,...,g, (3.2) 

M 

2g+i-] 

drj = ^ {k + 1- j)ak+i+j—dX, j^l,...,g. (3.3) 
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Figure 1: A homology basis on a Ricmann surface of the hyperelhptic curve of 
genus g with real branching points Ei, . . . , -E2g+2 = oo (upper sheet). The cuts 
are drawn from £^2i-i to E2i for i = 1, . . . ,g + 1. The b-cycles are completed 
on the lower sheet (the picture on lower sheet is just flipped horizontally). 



We denote periods 



2uj = { (t dhi\ , 2oj' = id) dhi 

'ofc / i,fc=l,...,g yJbk / i,k=l,...,g 



(3.4) 

2r]={- j drA , 2r]' = i- i dn) 

V Jak / i,k=l,...,g V Jbk / i,fe=l,...,5 

The half-periods u),u',r],r]' satisfy the generalized Legendre relation 
where Ig is the unit g x g matrix and 971 is the 2g x 25-matrix 



"-9 

on- 



V V 

Denote r = a;~^a;'-matrix belonging to the Siegel upper half-space = {rlr-^ = 
T, Im r > 0} and necessarily symmetric matrix k ~ rj(2Lu)~^. Jacobi variety of 
the curve is defines as the factor over the normalized period lattice, Jac(X) = 
C/lg ® r. Abel map 2t : {P e X} — > Jac(X) maps the point P of the curve to 
Jacobian by the rule 

dv, 

'Po 

where Pq is the base point of the Abel map. It maps a positive divisor D of 
degree n, D = Pi + . . . + Pn as follows 



2l(P) = 

Jp„ 



2l(£)) = 2l(Pi) + ... + 2l(P„). 
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In what follows we deal with non-special divisors (we remind that a special 
divisor D is a divisor for which it exists a meromorphic function on X with poles 
at most in D). The base point of the Abel map will be taken as Pq — (cO; oo) 
unless a contrary is stated. 



3.2 ^-functions and cr-functions 

Introduce canonic 0-function, 6{z]t), by the formula 

9{z;t) ~ exp {mrn^Tm + . (3.6) 

Fundamental Riemann vanishing theorem says that if for a vector z € ,Jac{X) 
one has that 9{z;t) ^ 0, then there exists a non-special divisor D of degree g 
and a vector Koc (vector of Riemann constant with the base point Pq = (oo, oo)) 
such that 

z^^2^{D)-Koo■ 

Moreover Riemann 0-function 

eiP)=0imiP)-miD)+K^;T) 

vanishes precisely in g points of the divisor D. Let (Ai, /ii), . . . , (Ag, /ig) be a 
non-special divisor. Denote 

Then the cr-function is defined by the formula 



det{2uj) 

X e{{2Ljy^h\T)exp{h' r]{2u;y-h;T). (3.8) 



l<i<j<2g+l 



The function (|3.8|l represents a g— dimensional generalization of the elliptic cr- 
function. It is invariant with respect to the action of certain subgroups of the 
symplectic group. It also has the following transformation property with respect 
to shift of the argument by a period, 

a{h + 2ujn + 2uj'm) = a{h)exp {{2rin + 2r]'m)'^{h + ujn + u)'m)] , (3.9) 

where n, m are integer t;- vectors. 

Variety of zeros of the cr-function called cr-divisor Q, coincides with the 
0-divisor (9). The cr-divisor in rational limit is given as zeros of a Schur- 
Weierstrass polynomial analyzed in 14 . 

Denote Abelian image 8i of the curve X, 



Qi^ Ive Jac(X)| j 



dh, Pex 
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At genera g > 1 the variety ©i belongs to the cr-divisor, Qi C 6. We shall need 
the following result from |37| describing vanishing of certain partial derivatives 
of the (T-function. 

Lemma 3.1 (Onishi, 2004) Let aj be partial derivative 

je3 ^ 

where 3 be multi-index given for different genera g of hyperelliptic curve X in 
the Table 1 
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2 


3 


4 


5 


6 


7 


8 


9 




3 




2 


2 


24 


24 


246 


246 


2468 


2468 





Table 1. Multi-indices corresponding to different genera g 



Then in the vicinity of the point h = (0, ...,0)^ G Qi the following ex- 
pansions are valid for the X-coordinate of the curve X and partial derivative 

A = 1 + 0(1), (3.10) 

CT3(/l)2 0(^29+2). (3.11) 

Moreover all derivatives crg(/i) vanish on /i G ©i for all multi-indices 3 C U, 
131 <P|. 

The proof of the lemma is based on the detailed analysis of the Schur polyno- 
mials associated with hyperelliptic curves and certain results of |14| . 



3.3 Kleinian p and (^-functions 

Introduce dimensional ( and p-functions as logarithmic derivatives of the 
cr-function as follows 

d 

Ci{h) = —\oga{h), i^l...,g, 
oui 

92 

Pt,]{h) g^_g^, ^ogcr{h), i,i = l,...,g, 
Pi,j,k{h) = ^ g^.g^ .g^^ ^og a{h), i, j, k ^ 1, . . . , g, etc. 

These functions have the following periodicity properties. Let n,m E Z,^ two 
arbitrary integer vectors. Then 

Cfe(M + 2a;n + 2a;'m) = Cfe(M) + 27;n + 277'm, k^l,...,g, (3.12) 
p3{u -\- 2ujn 2uj'm) ~ p3{u), (3.13) 
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where 3 be arbitrary set of two, three, etc., indices. 

Consider the principal object of the hypcrelhptic theory the master polyno- 
mial 

a 

F(A;Ai,...,A<,)- J](A-Afe), 

k=l 

where Ai, . . . , Ag are first coordinated of non-special divisor of degree g, Pi = 
(Ai, /xi), . . . ,Pg = (Ag, /ig). O.Bolza Jl^ found expression for the master polyno- 
mial in terms of p- functions (se also EJjEIjCSI for derivation and applications) 
in the form 

F{\- Al, . . . , Ag) = ^(A, h) = \^- Pg,g(/l)AS-l - ... - Pg^l(/l). (3.14) 

In what follows we shall refer to this parametrization of the master polynomial 
as to Bolza polynomial. 

Zeros of Bolza polynomial J'(A, h) are known to solve the Jacobi inversion 
problem, 

Al + . . . + Ag pg,g(/l), 

: (3.15) 

Al...Ag = (-l)3 + Vl,9(/i), 

whilst the /x-coordinates of the divisor Pi — (^i, Ai ),..., Pg = {fig,Xg) are given 

by 

d 



dug 



, fc-l,...,<?. (3.16) 



Among various relations of the theory we shall use, for our derivation the 
following formula for ^-function, 

B2g + 1 \ CXD / 

The derivation of this formula is documented in T and |12l IT^ and based on 
the finding of explicit expression for the fundamental bi-differential - so called 
Bergman kernel in modern terminology. 

Another set of formulae describe Jacobi variety Jac(X) as algebraic variety 

. g(g+l) 

in C^"*" 2 . The functions pggi and pit are related by 

PggiPggk = ^^PggPgiPgk - 2(pgipg_i,fc -|- pg,fepg-i,i) 

+ 4(pgfepg,i_l -I- pgipg^k-^l) + 4:pk-l,i-l - '^{pkA-2 + Pi,fc-2) (3.18) 

+ a.2gpgkPgi H "^l^i^igPkg + ^kgPig) + 

where ^ 

Ci,k = a2z-2Sik + -(a2i-i4,i+i + a2k~iSi.k+i)- (3.19) 
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We shall also use derivative of these formulae, 



Pgggi = (6pgg + Ol2g)pgi + Qpg,i~l " "^Pg-l,! + ^(5g,jQ!2ff-l, (3.20) 

which describe according to the KdV hierarchy. We remark that the above 
formulae were derived for the case ultraelliptic curves g — 2 m ^ where re- 
lation to integrable equations like KdV or sine-Gordon had not been noticed. 
The comprehensive generalization to hyperelliptic curves of higher genera was 
obtained to the first time in |12[ [T5| . 

We complete our brief description of u-functions by the addition theorem 
point+divisor for hyperelliptic cr-functions. 



Theorem 3.2 Let 

n(A,M) 



dh e Oi 



oo,oo) 



be Abelian image of the curve X and let h be a vector in general position in 
Ja,c{X). Then the following formula is valid 

where J'(A, h) is the Bolza polynomial and the multi-index 3 is given in 

the Table of Lemma \3.1\ 

To the best knowledge of the authors the formula H3.21|) as well the foregoing 
relation (|3.22() are new. We shall publish detailed proof elsewhere and give here 
only its short version. This proof is based on the Riemann vanishing theorem 
and solution of the Jacobi inversion problem in the form 13.15|1 what leads to 
the Bolza polynomial in the right hand side. Lemma |3. II permits to check the 
leading terms of the expansions in v near the point v = (0, . . . ,0)"^ G 8i in 
both sides of the equality. 

Corollary 3.3 Let fl ~ Sfcei^fc- where \I\ ~ g be non-singular even half- 
period represented in the form = 2ujn + 2Lu'n' with n,n' G . Then the 
following formula is valid 

''^'"^^f ^Oo(A)TT(A-^0, t;eei, OeJac(X), (3.22) 

where %{\) = (T(f2)2exp{A^ d^} with A = 2?7n -t- 2r^'n' , E, are branch 
points of the curve X and 3 is the set of multi-indices given in the Lemma \3.1\ 



4 (j-functional realization of the basic functions 

Consider /i-vectors of the form 

h = ixeg -n, eg = (0, . . . , 0, 1)^, (4.1) 
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where fl = (fii, . . . , ^g)^ be non pure imaginary even non-singular half-period 
supported by g branch points 

fl = ^2t„ |I|=.9. (4.2) 

Then the Bolza polynomial is identified as the S{x, A) function introduced in 
Section 2, 

T(A;/i) = 5(x,A). (4.3) 



4.1 Its-Mat veev theorem 

The following theorem describes algebro-geometric solutions of the Schrodinger 
equation with finite-gap potentials and was proved by Its and Matveev in 1975 
|27|. Different variants of proof the Its-Matveev theorem are known see e.g. [7], 

ca. m 

Theorem 4.1 (Its-Matveev, 1975) Let X be non- degenerate hyperelliptic curve 
of genus g. Let fl = (fii, . . . ,^g)^ be non pure imaginary even non-singular 
half-period supported by g branch points. Then the smooth and real potential 
u{x) and one-valued on the curve X Bloch function iplx^P), are given by the 
formulae 

U{X) ^ PggiiXBg -^n), (4.4) 



,Pix,P)^C{P) 



woo 



IXBn 



n 



xexp < la; 



dr„ - 



dr 



(S2g+1.0) 



(4.5) 



where 3 is the set of indices given in the Table of the Lemma \c/.U C{P) is a 
normalization constant which depends on the point P G X , eg = (0, . . . , 0, 1)"^. 



4.2 (T-functional description of the Weyl function and the 
wave number 

Theorem 4.2 Let conditions of the Theorem \4. l\ are satisfied. Then Weyl func- 
tion w{x, P) and wave number k{P) are given by the formulae 

W{X, P) - Cg '^XEg J dh\- CgiiXBg -ft)- ^ dr g . (4.6) 
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The wave number or quasi- momentum is 



k{P) = Ug UxBg - n + j - CgiiXBg ~ ^ ^ 

^ V 00/ (B2g+1,0) 



(4.7) 



Remark that due to the formula (|3.17() and mentioned identification S{x,\) as 
the Bolza polynomial the quasi-momentum is given as the second kind Abelian 
integral H2.12|l . Denote a-periods of k{P) as 



3C, = * d3C, i^l,...,g. (4.8) 



4.3 Normalization of the Bloch function 

Different normalizations of the Bloch function are accepted in various problems 
of physics. In our case the normalizing constant C{P) is computed from the 
condition (|2.14|l . We shall do that on the basis of the relation H3.21|l . As the 
result we obtain the following expression for the normalizing constant 

^^^^ " ^{s{x,\)y ^^-^^ 

The proof in general case is based on the addition theorem 17. Expression 
H4.9|l was alredy given in [7|, Chapter 8. We emphasize that it is valid both for 
ergodic and periodic cases. 

Further we shall fix the value (|4.9() for the normalization constant and con- 
sider only normalized Bloch functions ip{x^ P). 



5 Periodic and ergodic finite-gap potentials 

We shall now consider two classes of finite-gap potentials: periodic and ergodic. 



Definition 5.1 We shall call the vector U = {Ui, . . . , UgY" , 

U={2uj')-^eg (5.1) 

as the winding vector 

Definition 5.2 The finite-gap potential is called periodic if the curve X admits 
existence of an integer vector n G and real number U G R ?7 7^ such that 

U = nU, (5.2) 

otherwise the finite-gap potential is called ergodic. 

We remark that in generic case the finite-gap potential is ergodic, it is periodic 
only under special conditions on the curve presented above. 
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5.1 Condition of periodicity of the finite-gap potential 



The direct consequence of the above definition is that the finite-gap potential 
is periodic if and only if period matrix luj' satisfy to the equation for certain 
vector n G and real number [/ 7^ 0, 



1 



(5.3) 



One can easily see that under condition (|5.3|l the functions in x gki{x) 
pki i^xEg — Jl) are periodic in x with period i/U. Indeed, 



pki [i[x + 



U 



eg-fl] = pki (ixBg - H - 2uj'n) = pki {ixBg - O) 



Conditions of periodicity 1)5. 3(1 represent g transcendental conditions on the 
moduli of the curve. But these conditions are weaker then conditions of double 
periodicity in x of the finite-gap potentials. The last conditions are formulated 
in terms of period matrix of the curve r as follows 

Theorem 5.1 (|7|) The finite gap potential u{x) is elliptic function of x if and 
only if 



there exist such homology basis S = (oi, . 
matrix t has the form 



( Til k/N 

k/N * * 

:)< :)< 



, a^; bi, . . . , bg) that the period 
... 0\ 



(5.4) 



y :)< * >|t 

withk€{l,...,N-l}, iV>l,iVeN. 
• in the homology basis H the winding vector U is of the form 

J7= (*,0,...,0)^. 



(5.5) 



Conditions H5.4|l . (|5.5(l represent 2g — 2 equations and therefore the associ- 
ated double periodic potentials can be included as particular cases into the 
set of periodic potentials. These more particular potentials permit explicit 
analytic description and we shall use this circumstance in what follows. We 
remark that the reduction theory of Abelian functions to elliptic functions goes 
back to K.Weierstrass and A.Poincare and is exposed in the classical A.Krazer 
monograpghj3D|; the modern exposition and applications to integrable system 
are given, in particular, in 111,0; the Theorem 15.11 was recently considered in 
|22|. where an alternative proof is given. 
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5.2 Calculation of averages 

The average 



(5.6) 



is sensible to commensurability or non-commensurability of the frequencies 
(components of the winding vector U). In what follows we shall specify co- 
efficients of the polynomial (|5.6() as s^^"* or sj'^'' with superscript (p) taken for 
periodic case and (e) - for ergodic, we shall also provide {S{x, A)) with sub- 
scripts, {S{x,X))p and {S{x,\))e- 

Theorem 5.2 (Averaging in ergodic case) Let the components of the wind- 
ing vector U be all non- commensurable. Then the average {S{x, A))e is ergodic 
and given by 



{S{x,X)), = XS + J2^' 



-1 Je) 



where 



det2w' 



-det 



In particular, for g = 1 



{S{x,\))e=X + 



(5.7) 



and for g — 2 



{S{x,X))e 



det 



A^ + A- 



det 



^11 


^12 


V21 


V22 




^'12 


^21 


^22 



det 



det 



V21 


V22 


^21 


^22 


^'11 


^12 


^21 


^'22 



In the case of commensurability we obtain another result for the average. 

Theorem 5.3 (Averaging in periodic case) Let the components of the wind- 
ing vector U be all commensurable. Then the average {S{x, X))p is periodic and 
given by 



{S{x,X))p = zS + 2_^z- 



(5.8) 
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wher 



( <,i 



1 

det2w' 



-det 



'3,1 



1. \ 
J 



"g-1,9 



(5.9) 



In particular, for g — I 

and for g — 2 

(5(x,A))p = A2 + A- 



I ,y 



det 



^11 ^12 

%1 %2 



det 



det 



11 



■^12 



det 



^^21 


^22 


<1 


'?12 


^'11 


^^12 


^'21 


^^22 



We see that for the case of elhptic curve the averages in ergodic and periodic 
case coincide, but for genera bigger than one {S{x, X))p and {S{x, A))e are given 
by different expressions. The comprehensive proof of these formulae wiU be 
given elsewhere. 

6 Wannier function 

Definition 6.1 The Wannier function Wn{x) in the n-th spectral band, n = 
1, . . . , g, is defined as the following integral of normalized Bloch function, 



1 



''^ n J a 



■0(a;, A)d3C(A), n=l,...,g, 



(6.1) 



where 3C„ are an-periods of the differential of quasi-momentum \2.1!^) . 



Using translation operator one then construct a countable set of Wannier 
functions, 



with 



1 



V'i(a;, A)dDC(A), n = 1, . . . , g, Z e Z 



(6.2) 



I CT (^/^ d/i - ixBg + + 2in' 

V^SO^iAjy crj (^/^ dh^ a{ixeg il- 2m') 

(p p 
ix J drg ~ {2m' + flf J dr 



(6.3) 



(S2g + 1,0) 



(£^29+1.0) 
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where fl' is a pure imaginary non-singular even half-period. In the case of genus 
g — 1 the above definitions are those given in |10| . 

To complete this definition we must compute periods H2.13|l . To this regard 
the following lemma apply. 

Lemma 6.1 Denote and Sk are coefficients of the polynomial 
{Six,X)) =\^ + Sg\a-^ + ... + si. 

Then 

f ^1 \ ( 2xg,i -si. ^^ 

(6.4) 







1 2Xg^l - Si 


V ) 







In particular in the case g = 1, si = rj' joj' , >t\\ — >t — rj/2tjO and 



ITT 



3C1 - — . (6.5) 

Note that in the case g = 1, the Wannier function H6.1|l coincides with the 
classical one derived in [TDj . 



/ y \ 1/2 r-Tv/T 



The principal expected properties for the set of Wannier functions, Wk{x)^ 
fc = 1 , . . . , g is orthogonality 

w''P{x)wP{x)dx = 5ki6k'v, k,l^l,...,g, k'J'eZ (6.6) 

and completeness in space of eigenfunctions of the Schrodinger operator with 
finite-gap quasi-periodic potential. We shall postpone the discussion of these 
questions to a future publication. 

Since the normalization of Bloch function is different for periodic and ergodic 
cases we shall distinguish two kind of Wannier functions - the classical Wan- 
nier functions associated with periodic potentials and the quasi-periodic Wan- 
nier functions associated with quasi-periodic potentials. Both kinds of Wannier 
functions are defined by the same formula involving Bolza polynomial in the 
expression of the normalization constant. Normalizing constants, however, are 
computed in different ways for periodic and ergodic cases. 

6.1 Power series for the Wannier function at |x| ~ 

Theorem 6.2 The Wannier function of the first lower energy band Wi{x) ad- 
mits the following expansion 

Wi{x)-f:i^W^"'^^''- (6-7) 
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Here 

W^f^^ =^Mfegpfe, (6.8) 
fe=i 

the quantities Mk are given as a-cycles of second kind Abelian integrals, 



1 / xfc M'' + ^9^'^^ + • • • + Si 



M, = ^f \- ^ oT^^^' (6-9) 

T/ie quantities qti are coefficients of polynomials in z, Qp — X]f=o defined 
by recurrence 

m=0 ^ 

wzi/i conditions 



2p + 2 



r/ie /irst /ew coefficients are 
VFf ) = Mo, 



W^^'> = M3 + 6pg3(n)Af2 + (l2Pss(fi)2 + 14pgggg(rj))Ml, 
+ (2p333333(f^) + 28pgggg{n)Pgg{n) + 8 Pgg ( ) 3 ) Mq • 



Here in its turn 

Pggggin) = (6 ^ i;, + ftsfc ) ^ + 4 ^ i;, i;^ , 

iei iei ijfcei 

where we used H3.2()(l . 

We remark that the proof of this theorem requires an expansion of the ex- 
ponential in X which is a generating function for Schur polynomials '3"2! and 
subsequent integration of these polynomials (details will be given elsewhere). 
The integrand can be expressed recursively in terms of second kind Abelian dif- 
ferentials described below. In the case of genus one curve and elliptic periodic 
potentials these integrals are given as standard hypergeometric functions [IT)] . 

Another remark is that in the course of the proof a reciprocal hyperelliptic 
curve X"^ of the same genus and dual to the initial hyperelliptic curve X, 
naturally arises. 
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Definition 6.2 We shall call reciprocal curves X* , k = 1, . . . ,g associated to 
the curve X the set of hyperelliptic curves of genus g given by the formulae 

/x#'=4(Af + .s<,Af-i + ... + si) H (X-E,) 

j:e{i,...,2ff+i}-ifc (6.10) 

= 4(A2^+i+a|A2^ + ... + «„#), 

where Ik is the set of g indices, Ifc C {1, . . . , 2g + 1}. 

We remark that the ampHtudes of the Wannier pulses W„(0), n — 1, . . . , g 
are given as periods of meromorphic hyperelliptic integral of corresponding re- 
ciprocal curves. 

To describe canonical differentials on xf it will be convenient to introduce 
the set of second kind differentials, which we shall call the Wannier differentials. 



dWi^=A"J= — dA, n = 0,.... (6.11) 

V ll»e{i,....2g+i}-i(.-^ ~ ^i) 

Canonical holomorphic differentials dh* = {dhf , . . . ,dh*)'^ of the curve are 
given as 

d/if = ^dW(") = ^dA, fc-l,...,<?. (6.12) 

The fact that holomorphic differentials of the reciprocal curve are gives as deriva- 
tives in moduli of meromorphic differentials resembles principle feature of the 
algebraic curves appearing in Seiberg-Witten theory j^jjlUl, see also [23 where 
various forms of Seiberg-Witten curves were discussed. 

Introduce also canonical second kind differentials dr* — {drf , . . . , drf)"^ as 
in (E31) 

29+l-i ^fc 

drf - E (fc + l-J>f+i+,^dA j = l,...,5. (6.13) 

k=j ^ 

The first g differentials dW^*^^ , k ~ 0, . . . , g — 1 are expressible linearly in terms 
of 2g differentials dr* and dh* . Periods of the differentials dW^"-^ with higher 
n are also expressible in terms of aforementioned 2g differentials due to the de 
Rham theorem |25| . 

The last point which we wish to discus here is the asymptotic expansion 
of the Wannier function at x — > -|-oo. To evaluate this expansion one can 
use a variant of the steepest descents method applicable to the case when the 
integrand contains a multiplier vanishing in the saddle point. The following 
theorem apply. 

Theorem 6.3 At x —i- oo the Wannier function of the lower energy band 
[Ei,E2] for the one gap potential has the following asymptotic expression 



W^{x) Re I J MM / J](Ao - £;,)V(x, Ao) 



2(5'(x,Ao)) 



(6.14) 
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where Ao is a solution of the equation 

{S{x,X)}=0, (6.15) 

and Qo(A) is given in the Lemma \3.!A 

The formula for asymptotic expansion naturally generalizes the result for 
genus one derived in ^Ul and shows the universal character of the decreasing 
low at infinity of the Wannier function: exp {cqx} x"'^/^ with phase cqx in the 
saddle point. 

7 Conclusion 

In this paper we have defined Wannier functions associated with permitted 
zones of quasi-periodic finite gap potentials. A number of interesting questions 
remained beyond our considerations. One question is whether Wannier function 
introduced above generate an orthogonal and complete set of in the space of 
eigenfunctions of the Schrodinger operator. Moreover, a discussion of the first 
non-trivial examples, as well as a comparison between our analytical expressions 
for amplitudes and asymptotic behaviors of Wannier functions with numerical 
calculations as done done in |E3 fo'" the one gap case, would also be desirable. 
We plan to present this material and other related problems in a forthcoming 
detailed publication. 
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